q- Functional Field Theory for particles with exotic 

statistics 



K.N.Ilinski x ' 2 ' 3 * G.V.Kalinin 1 t and A.S.Stepanenko x ' 2 ' 4 * 

1 IPhys Group, CAPE, St. -Petersburg State University, 
14-th line of V.I. 29, St- Petersburg, 199178, Russian Federation 

2 School of Physics and Space Research, University of Birmingham, 
Edgbaston B15 2TT, Birmingham, United Kingdom 

3 Theoretical Department, Institute of Spectroscopy, Troitsk, 
Moscow region, 142092, Russian Federation 

4 Theoretical Department, St. -Petersburg Nuclear Physics Institute, 
Gatchina, St. -Petersburg, 188350, Russian Federation 

Abstract 

In the paper we give consecutive description of functional methods of quantum 
field theory for systems of interacting (/-particles. These particles obey exotic statistics 
and appear in many problems of condensed matter physics, magnetism and quantum 
optics. Motivated by the general ideas of standard field theory we derive formulae 
in g-functional derivatives for the partition function and Green's functions generating 
functional for systems of exotic particles. This leads to a corresponding perturbation 
series and a diagram technique. Results are illustrated by a consideration of an one- 
dimensional (/-particle system and compared with some exact expressions obtained 
earlier. 



1 Introduction 

In this paper we deal with so-called g-particles, i.e. particles which appear as a result of 
quantization of a Hamiltonian classical dynamics on g-deformed graded-commutative alge- 
bras [1]. Their creation and annihilation operators a\,ak obey the following commutation 
relations: 

OfcOj = KqkjCijak , a k Oj = Kqkj0^a\. , a^a] = Kq^a^ak + 5kj , 

ctkCLk = a\a\ = for g-fermions . (1) 
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Here the deformation parameters qkj possess the property 

qkj = e^ fcj , 4> k j = -4>jk , 4>kj e JR 

and n serves to unify formulae for deformed bosonic and deformed fermionic cases. As usual 
it has the form: k = +1 for g-bosons and k — — 1 for g-fermions. 

The notion of the exotic quantum statistics seems to be rather artificial or mathematical. 
Number of papers where particles with the statistics emerged in a context of the parastatis- 
tics [1, 2, 3], g-extended Supersymmetry [4, 5], Parasupersymmetry [4, 6] and other similar 
problems apparently prove the statement. However, solid state physics, quantum optics 
and theory of magnetics give examples of other kind. Indeed, anyons (particles with exotic 
braiding statistics) are important in some attempts to understand the physical features of 
planar systems [7], the quantum Hall effect [8] and high temperature superconductivity [9]. 
In contrast to these examples where anyons serve as auxiliary objects for the construction 
of one of possible scenarios there is a wide field in the quantum nonlinear optics where 
g-particles are the main components. This is a theory of the collective behavior of exci- 
tons with small radius (Frenkel Excitons and Charge- Transfer Excitons (CTE)) [10]. The 
studies investigate possibilities of formation of the Frenkel biexcitons and the observation 
of phase transitions in exciton systems in molecular crystals (Bose- Einstein condensation 
of excitons [11], structural phase transitions in crystals with high excitonic concentrations, 
dielectric-metal phase transition in a system of interacting CTE [12] and others). Strictly 
speaking excitons are not particles. They are quasiparticles describing molecular excitations 
and are of great importance in the analysis of nonlinear optical processes which accompany 
propagation of high-intensity light fluxes whose frequencies are in the range of the exciton 
absorption bands [13]. These excitons obey exotic statistics (Pauli statistics) [14] coinciding 
with g-particles statistics for q = —1. The general case of q = e 1 ^ arises if we try to take into 
account phenomenologically some nonlinear effects (such as the difference in the creation 
time of molecular excitations for different types of molecules). This effect can be modeled by 
the change of the Paulion commutation relations to those of the g-particles using the method 
developed in [15]. Noteworthy, even the investigation of the behavior of low dimensional ex- 
citon systems is meaningful. The best example is the exact solution for one-dimensional 
Paulion chain [16] caused great advances in the theory of the so called J-aggregates, i.e. 
molecular aggregates with unusually sharp absorption band ([17] and Refs. therein). The 
investigations of exciton systems on interfaces closely connect with the successes of contem- 
porary technology. All these show that g-particles find deep applications in modern physical 
theories and motivate our objective to derive the appropriate field theoretical technique for 
them. 

Recently, it was shown that g-functional form of Wick's theorems for creation and an- 
nihilation operators of the g-particles can be formulated and they have the same formal 
expressions as fermionic and bosonic ones but differ by a nature of fields[18]. It means that 
in the case of the g-particles certain g-deformed algebras should be used exactly as it was 
with Grassmann algebra in the case of fermions or the complex numbers in the case of 
bosons. This fact allows us in the present paper construct a machinery of the quantum field 
theory for the exotic particles going along the way of standard textbooks. 

In a sense, the present work may be considered as a consequential step to the quantum 
field theory for the exotic particles which follows from previous papers [1, 19, 18]. Indeed, 
a construction of classical and quantum dynamics on graded-commutative spaces [1] stated 
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a connection of the g-particles and the g-deformed classical variables. These variables then 
have been used to introduce the corresponding coherent states and to derive (/-functional 
integrals [19] for g-particle systems. It is well-known that Wick's theorems present other 
way to the functional integrals and field-theoretical methods. Since the g-functional Wick's 
theorems have been proved [18] we now make the ends met and derive the field-theoretical 
technique and functional integrals from the theorems. This step (this paper) completes the 
program. 

The paper is organized as follows. In the next section we give a functional represen- 
tation for the partition function and Green's functions generating functional. Then in the 
sections 3 we illustrate the developed technique with an example of calculations for simple 
a one-dimensional g-fermionic system. Earlier, in Ref. [20] we presented exact results for the 
partition function and Green's functions for the system using a functional integral method. 
This gives us a possibility to compare these two approaches and to state their consistency. 
Closing remarks conclude the paper. 



2 (/-Functional field theory 

In the quantum statistics all equilibrium physical properties are described by the density 
matrix given by the operator 

p = exp[-f3H] , f3 = 1/kT . (2) 

In particular, the equilibrium thermodynamics is governed by the partition function which 
is a trace of the density matrix: 

Z = tip . (3) 
Mean value ((b)) for an arbitrary quantum operator b then may be calculated as 

(($)) = Z-hr [pb] . (4) 

In general, the partition function cannot be calculated exactly and one should develop a 
perturbation theory. To do this, as usual we need to divide the Hamiltonian into two 
parts: so-called "free" Hamiltonian H which may be treated exactly and an interaction V 
considered as a perturbation: 

H = H + V, ifo = £ e *4afc ■ (5) 

k 

Since there is no appropriate way to treat hopping terms exactly we include them in the 
interaction. 

The ratio of p-matrices Po 1 p = exp[/3i?o] exp[— (3H\ coincides with the Euclidean evolu- 
tion operator U (j3, 0) in the interaction representation. Hence, the partition function and 
various averages require explicit expression for the operator U (/3, 0) in some approximation 
and they may be written as 

Z = tr p = tr [p U(P, 0)] = Z ((U(P, 0))) , (6) 

where 

((b)) = Z-Hr [p b], Z = tr[p ]. (7) 
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Creation and annihilation operators a\a in the interaction representation take the form 



4(r) 



e THo ale- THo 



a k (r) = e THo a k e- THo 



(8) 



From (5) we get an explicit form for the evolution of the operators: 



a' k [T) = ale 



a k {r) = a k e 



So the problem is reduced to the calculation of the evolution operator U{ji, r 2 ) in the inter- 
action representation. It is defined by the following equation: 



(9) 



with the Hamiltonian in the interaction representation 

V(t) = e H » T Ve~ H ° T . 
The solution of eq. (9) with the initial condition U (r, r) = I is the Volterra series: 
U(t u t 2 ) = e^ ori e^ (r2 - ri V OT2 = 

= E(-l) n r<i*i... [ Tl dt n 9(l...n)V(t 1 )...V(t n ) . 

r. J To J TO 



n=0 



(10) 



If the operator V is an operator functional of the bosonic type then RHS of the equation 
can be represented as the standard Dyson To-exponent by the symmetrization with respect 
to the permutation of time variables t k . In contrast to the undeformed case this condition 
holds only for very special cases. In general, the symmetrization of RHS give us another 
type of T-exponent. But this type of T-exponent is not consistent with nature of operators 
which makes difficult to deal with it. It would be natural to try to g-symmetrize RHS of 
(10) but different monomials in the operator V are permuted in different ways. So we can 
not g-symmetrize the Volterra series with respect to the permutation of the whole operator 
functionals V (or time variables) and we do not obtain T^-exponent. But we can deal directly 
with Volterra series (10) because each term in the sum can be represented as T^-product 
(due to a presence of the ^-function) [18]. 

Let us now reduce the Euclidean evolution operator U(ti,t 2 ) to the normal form. We 
suppose that the operator V is g-symmetrical one, i.e. V =Sym q V [18] and it does not 
contain any time derivatives. Adding the sign of the g-chronological product to RHS of 
eq.(10) and applying Theorem 3 of Ref [18] we obtain the following rule of reducing of 
eq.(10) to the normal form: 





[5 A 5 ] 


exp 


da da + 



U(a + ,a;T 1 ,T 2 ) 



,.+=«t 



11^ 



where following [18] we have introduced classical variables a, a + corresponding to the oper- 
ators a, a) and satisfying the following permutation relations 



a k aj = K,q kj aja k , a^aj = Kq kj ajal , a k aj = nq kj aja k 
a k a k = a k a k = for g-fermions . 



(12) 
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In formula (11) the g-chronological contraction A is defined by the relation 



(13) 



In the form J^At^j- the summation over discrete variables and the integration over contin- 
uous ones are implied. 

Now we encounter the problem of calculating of the mean value of an operator standing 
in the normal form 

«NF(ot,a)» =F(£,^)((Nexp[at C + c + a])) . (14) 

We imply the auxiliary fields (sources) c + ,c obey permutation relations (12). 
So the problem is in calculating of the following object: 

/(c + ,c) = ((Nexp[a f c + c + a})) = ((exp(a f c) exp(c + a))) 

It is obvious that the problem falls into a set of one- dimensional ones: 

fk(at,a k ) = Z^hr [exp(-/?£ jfe aj.a jfc ) exp(ajc fc ) exp(c£a fc )] 
Z k = (1- Kexp{-(3e k ))- K 

The calculation reveals the results of the undeformed case and we get the following expres- 
sions: 

exp(-/3e k ) 



fk(4,c k ) = exp(Kn k ata k ) , n k = , 

1 — /texp(— (3e k ) 

Then for any operator functional F it is possible to write the following relation: 



((NF(a\a))) = Fi 



8 d 



ti- 



de ' dc + 



i exp[c + dc] 



(15) 



(16) 



c+=c=0 



where F(a) means that in each monomial of F multipliers stand in the inverse order and 
the following notation is introduced: 



Using the identity 



d d 

exp[c + <ic] exp[£ + c + c + £] = exp[— d— — ] exp[£ + c + c + £] 



(17) 



(18) 



we can rewrite (16) as follows 
«NF(at,o))> 



= ex 4^ ( l- K ^ )exprc+c+fl 



c+=c=0 



c+=c=0 



(19) 
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The second equality is due to the fact that d contains 5-symbol (it is permuted without any 
phase factor). Let us introduce S-matrix functional as 



R(a + , a) = exp 



U(a + ,a;f3,0) 



(20) 



At collecting together eqs.(6, 11, 19) we get final expression for the partition function: 

Z/Z = R(0) . (21) 

In the last formula the deformation parameter q appears only in permutation relation for 
variables (12) and derivatives. 

Now we consider an application of the Wick's theorems to the calculation of Green's func- 
tions generating function. S-matrix Green's functions (without vacuum loops) are defined 
by the following relation 



H n (x!, ...,x n ) = ((T D [<£ H (zi), • • • , (fH(x n )])) 



(22) 



where x = (k, s,t), (p means a) or a (ip(k, l,t) = a k (t), ip(k, 2,t) = d^it)) and the subscript 
"if" means that the operator is in the Euclidean Heisenberg representation 

n (x) = ip n (k, s, t) = e tH 0(k, s, 0)e~ tH . 

Using the group property of the evolution operator we obtain 

^h(xi) • . . . • 0hM = U-\(3, 0) U((3, fi)0(xi)l7(fi, f 2 ) . . . <p(x n )U(t n , 0) . 

It means that we can write down the following equality: 

{(T D [U(/3, t^x^UiUM) ■ ■ ■ 0(x n )U(t n , 0)])) 



((T D [^ H (ll),...^H(ln)])) 



(W,o)» 



(23) 



Now we shall reduce the expression in the angle brackets to the normal form and then apply 
eq.(19): 



G n (x u ...,x n ) = ({T D [U(P, hWxjUfa, h)... 0(x n )U(t n , 0)])) 



= exp 



ld_ d_ 

2 d V 9 d V 



U (</?; (3, tiMx^Ufo; ti, t 2 ) • • ■ (p(x n )U((p; t n , 0) 



(24) 



where matrix temperature propagator has the form 

9 = 



n(d + A) 
d + A 



These functions G n (x±, . . . ,x n ) might be obtained from the generating functional for the 
Green's functions: 



G(A) = exp 



ld_ d_ 

2 dip 9 dip 



U(<p,A;P,0) 



(25) 



ip=0 



where 

U(tp, A; (3, 0) = U((p; (3, 0) exp ipA . 
Indeed, the Green's functions G n (xi, . . . ,x n ) are determined by the relation: 



(26) 



G n (X\ , . . . , X n ) 



5, 



5 t 



8A(x x ) ' ' ' 6A(x n ) 



G{A) 



(27) 



where the subscript "chr" means a specific procedure of the differentiation: each monomial 
term should be chronologically ordered, the variable to be differentiated should be moved to 
the most left position relative to other variables with the same time and then canceled. 

The formulae in the above paragraph are sufficient to work out a field theoretical tech- 
nique to calculate Green's functions and use various tricks of a field-theoretical machinery. 
However, we stop at this stage for a moment to establish a useful connections between the 
generating functional (25) and the S-matrix functional (20). All we need to do this is the 
following easy looking formulae: 



exp(-y^)F | A j exp^A) = fIa+-^- 



(28) 



exp\A-^\F(<p) = F(A + ( p) . 



(29) 



In these formulae we have right derivatives instead of left ones but it is not a problem to 
change all formulae in a proper way. Indeed, due to the locality of A and d we get: 



R((fi) = exp 



2 dip^ dip 



U(<p;P,0) 



(30) 



G(A) = exp 



nd_ d_ 

2 d<f^ dip 



U(<p,A;P,0) 



ip=0 



For the "free" theory (V = 0) using (28) we obtain 



(31) 



G^\A) 



exp 



nd_ d_ 

2 dip^ dip 



exp ipA 



exp 



<p=0 



K A A 

2 



(32) 



By the same trick it is possible to rewrite the expression for U(ip, A; (3, 0) as 



Uip. A: L0) = V (K^:.i()), : . 1 . 



(33) 



This gives us the possibility to present the generating functional in the form: 



O(A) = I |^;/3,0)oxi) 



(34) 
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From (31) and (28) we obtain 



G(A) = exp 




U(<p;P,0) 



(35) 



</?=o 



As it follows from g-symmetry property of A, d and definition (30) the generating functional 
has a very simple connection with the S-matrix functional: 



G(A) = exp (^AgAj exp 



d_ 

dp 



R(<p) 



<p=0 



or, using (29), in more compact form: 



G(A) = exp (%AgA) R{gA) = G^(A)R(gA) . 



The inverse relation is also useful and looks as (A — > g 1 A) 

R{A) = exp (-^<T^) Gig' 1 A) . 



(36) 



(37) 



Let us emphasize now that formulae (30), (31), (36) and (37) are absolutely identical to 
the corresponding formulae of the standard theories with the only difference in the nature of 
the fields. This is the central point to derive a proper diagram technique. 

However, before turning the attention to the diagram rules we spend a minute to show 
a connection with the g-functional integral formalism developed in Ref.[19]. The basis for 
the bridge is eq.(34). To this end we note that the Gaussian exponent in RHS of (34) has 
to be expressed in a g-functional integral form and then action of the differential operator 

U (k-^;/3, 0^ is processed explicitly under the sign of the functional integral by the usual 

way. This gives the complete action in the exponent under the g-functional integral and, 
as a result, the g-functional integral representation for the Green's functional generating 
functional emerges. It is exactly the same expressions which was obtained in paper [19] for 
a situation of an additional internal (anyonic) gauge field. 

Let us draw outline of a diagram technique. From (36) and (37) it obvious that a knowl- 
edge of the S-matrix functional R((p) is equivalent to a knowledge of all Green's functions 
G n and vise versa. So we consider here the diagram technique for the S-matrix functional 
only. 

From the definition (20) we have the following perturbation theory series for S-matrix: 



R((p) = exp 



ld_ d_ 

2 dp 9 dp 



n=0 



/ V . . . f P dt n 0(1... n)V(ip(h)) . . . V(p(t n )) , (38) 
Jo Jo 



It is convenient to calculate the S-matrix functional (38) in terms of diagrams. Each mul- 
tiplier V(p(tk)) is represented by a vertex on line (all the vertices are ordered in time). 

Action of -Q^g-§^ corresponds to adding of a line g connecting a pair of vertices. The line is 
added by all possible ways as each derivative ^ acts on any multiplier V(ip). In particular, 
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two derivatives of the quadratic form may act on the same vertex V . Such lines are called 
tadpoles. 

The result of action of the differential operation on a n-th term in the sum in (38) can 
be represented as a sum of diagrams consisting of iV time ordered vertices with any number 
of added lines. The vertex with n attached lines is associated with the following expression 



V n (x 1 , . . . , x n ; ip) 



8 n V{^) 



S(p(x 1 ) . ..5(p(x n ) 



(39) 



The arguments x of the multipliers (39) are contracted with the corresponding arguments of 
lines g. Multiplier V(ip) = Vo(ip) is called generating vertex. If the interaction is polynomial 
in fields then finite number of multipliers (39) are non-zero only. 

For a generic term in the sum (38) the following representation is valid [18]: 



cxp 



^ dpi 9 dip k 



ik 



-l) n /V . . . f<\t n 6(1... n)y(^i(*i)) . . . V(<p n (t n )) 
Jo Jo 



(p 1 =...=tp n =<p 

(40) 

The diagonal terms of the quadratic form in the exponent correspond to adding of the 
tadpoles. They can be accounted by introducing of a reduced vertex 



KedM = exp 
and, hence, (40) can be rewritten as 



id_ d_ 

2 dip 9 dip 



V(<f(t)) , 



exp 



1 



d d 
g- 



y 

2 ^ d<pi d<£ k 



(-1)" /%..- f P dt n e(l...n)V nA (<p 1 (t 1 ))...V nA (<p n (t n )) 
Jo Jo 



<Pl = ---=<Pn=<P 

(41) 

The remaining terms in the differential operation add lines between different vertices. V(ip) 
represents Sym 5 -form of the interaction, V icd (ip) does its N-form [18]. 

Expression (41) is a generic term of the perturbation theory series and can be represented 
in graphic terms (diagrams). Due to the time ordering of vertices diagram rules and, in 
particular, procedure of calculation of symmetrical coefficients, differ from the standard 
ones but, however, straightforward now and can be easily adapted for symbolic computer 
calculations. 



3 Illustration: one-dimensional g-fermion system 

In the previous section we develop a general technique for a calculation the partition function 
and the Green's functions generating functional. This section is devoted to an application of 
the technique. As a field of the application we choose so-called cyclic g-XX-chain which was 
introduced and exactly solved in Ref. [20]. There it was shown that the partition function 
and the two-point correlation functions for the model can be calculated in explicit form so 
now we can use them to test our results. Let us remind the Hamiltonian of the g-XX-chain: 

M 

Ho = B <Krfim , 
m=l 
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M-1 




V(a\ a) = A (a^a m +i + a} m+l a m ) + A(a ] M ai + a\a M ) , 

m=l 

where the creation and annihilation operators a^, a m (m — 1, 2 . . . M) obey the commutation 
relations 

a k a,j + ga 3 -a fc = , a fc aj + q~ l a)a k = , g = e * 27r//W , 1 < fc < j < n , , . 

a k a[ + 4«fc = 1 , [a\] 2 = [a k f = . 

The relations coincide with (1) in the case <pj k = 2nl/N for 1 < j < k < M. The explicit 
result for the partition function of the model has the following simple and compact form: 

1 + <? exp{-/?[5 + 2A cos(^(A; + ^^))] }) (43) 

We compare this formula with eq.(21) applied to the case of cyclic g-XX-chain. More precise, 
we compare the perturbation series in the hopping parameter A. 

There are several technical simplifications to note. First of all, the exponents e £fc ^ 2 ~* 1 - ) 
may be excluded from (13, 17) since a generic term Ilfc=i V{tk) contains an equal number of 
a + {t k ) and a{t k ) and the energy on all sites are equal to B. This leads to the cancelation 
of the exponents. Second point is the fact that for this particular case the deformation 
parameter q appears in formulae starting with the order M with respect to the interaction 
(hopping). Indeed, for lower orders the only terms contribute to the answer are the monomial 
terms with equal numbers of products a k ai and a^a k (k = I ± 1). When we collect such 
pairs of products together using the permutation relations no phase factor appears. This is 
because of the property of this products to commute with terms as a+a m (m ^ k,l). When 
derivatives from expression (20) act on fl a t a i a t a k the deformation parameter does not 
appear. The result is natural: m-th order corresponds to m hoppings. When m < M there 
is no cyclic path for a particle to go and the deformed statistics does not play a role [20] . 

Summing up, in the lower M — 1 orders the contributions are pure fermionic. It is 
interesting to compare nontrivial orders. For any particular finite M it is possible to do by 
a straightforward calculation. Here, for sake of simplicity, we present the simplest nontrivial 
case M = 3 since the case illustrate the situation by a very clear way. First nontrivial order 
for the example is third order in the hopping constant and it might be calculated immediately 
from eq.(21): 

-J- = 1 + 3(3 2 A 2 [n - n 2 } + (3 3 A 3 [-n + n 2 (2 + (q + q- l )/2) - n 3 (l + (q + q' 1 )^)} . (44) 

n is defined by eq.(15) and in this case is equal to e~ l3B /(l + e _/3B ). 

It is not difficult to see that the exact formula (43) gives the same result. Let us note 
that the expressions in the round parentheses [(2 + (g + q~ l )/2) and (1 + (g + g _1 )/2)] are q- 
symmetrical coefficients and play the same role here as the standard symmetrical coefficients 
do for the case q — 1. Technically, the parameter q appears here from the permutation of 
ordered vertices. 

4 Conclusion 

In the paper we derived the representation (in g-functional derivatives) for the partition 
function and the Green's functions generating functional starting with the Wick's theorems 
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for creation and annihilation operators of g-particles. Derived representations may be rewrit- 
ten making use of functional integrals over functions on (g-deformed) graded-commutative 
space. It is not difficult to see that this would bring us back to formulae of Ref. [19]. Diagram 
technique of the paper follows from both formalisms (functional derivatives and functional 
integrals) equally easy. All said above let us state that the quantum field theory machinery 
for particles with deformed exchange statistics is constructed. It is similar to the standard 
bosonic and fermionic theories and differs by the g-deformed nature of fields (i.e. correspond- 
ing classical variables) only. We checked the technique on an one-dimensional example. The 
next step (which is much more physically intriguing) is to examine a mult i- dimensional case. 
Long-range order, instabilities, influence of disorder might be some of questions of interest 
for a such investigation. We return to them in a forthcoming paper. 
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